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The weak radiative corrections to the Drell-Yan process above the Z-peak have been studied. 
The compact asymptotic expression for the two heavy boson exchange - one of the significant 
contributions to the investigated process - has been obtained, the results expand in the powers of 
the Sudakov electroweak logarithms. At the quark level we compare the weak radiative corrections 
£ — ' to the total cross section and forward-backward asymmetry with the existing results and achieve 

a rather good coincidence at y/s > 0.5 TeV. The numerical analysis has been performed in the 
high energy region corresponding to the future experiments at the CERN Large Hadron Collider 
£NJ ■ (LHC). To simulate the detector acceptance we used the standard CMS detector cuts. It was shown 

» I ' that double Sudakov logarithms of the WW boxes are the dominant contributions in hadronic cross 

f"S |. section. The considered radiative corrections are significant at high dilepton mass M and change 

■ the dilepton mass distribution up to ~ +3(— 12)% at the LHC energies and M — l(5)TeV. 

: 

(N ' I. INTRODUCTION 

> : 

Despite the fact that the Standard Model (SM) for more than twenty years has been considered a consistent and 
£f~1 • experimentally confirmed theory, the search of New Physics (NP) manifestations has still been continued. As the 
possible (and, at least, most popular) phenomena of NP, we can name the supersymmetry, the large extra dimensions 
[lj and extra neutral gauge bosons [2J. The light on this paramount problem of modern physics can be shed by the 
coming in the near future experiments at the collider LHC. 

The experimental investigation of the continuum for the Drell-Yan production of a dilepton pair, i.e. data on the 
f**^ 1 cross section and the forward-backward asymmetry of the reaction 

^. pp->7,z->z+rx (i) 

at large invariant mass of a dilepton pair (see [3] and references therein) is considered to be one of the powerful tool in 
the experiments at the LHC from the NP exploration standpoint. The cause of this is high sensitivity of the dilepton 
continuum to any modification of SM induced by NP. The research of the NP effects is impossible without the exact 
knowledge of the SM predictions including higher-order QCD and electroweak radiative corrections. 

There have been numerous publications on this problem; let us now determine the place of the presented paper 
(further we will discuss only the electroweak corrections, the QCD contributions are out of the presented paper). So, 
all of the electroweak corrections to ([1]) can be divided conventionally into three categories: 

• (I) electroweak corrections induced by gauge Boson Self Energies (BSE), 

• (II) the other QED corrections (i.e. radiative corrections induced by at least one additional photon: virtual or 
real), 

• (III) the other Weak Radiative Corrections (WRC) (i.e. radiative corrections induced by additional heavy 
bosons: Z or W), 

The first and the second contributions have already been studied (see papers on pure QED corrections and the 
QED and electroweak corrections in the Z-peak region and above in and numerous papers cited there). Naturally 
the contribution (II) also requires the consideration of the diagrams with real bremsstrahlung (to cancel the infrared 
divergence) and, therefore, causes the experimentally-motivated kinematical cuts for observed photon; all of that has 
been taken into account in paper 
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FIG. 1: Feynman graphs for the Born (a) and one-loop (b-h) diagrams with additional virtual heavy bosons corresponding to 
the WRC. Unsigned helix lines mean 7 or Z. 

To describe all contributions to electroweak corrections, it is convenient to use the terminology of the so-called 
electroweak Sudakov logarithms [a], i.e. the expressions which are growing with the scale of energy, and thus giving 
one of the main effects in the region of large invariant dilepton mass. By now extensive studies have been done in 
this area. For instance, the weak Sudakov expansion for general four-fermion processes has been studied in detail 
(see, for example, 0] and the recent paper Q along with the extensive list of references therein). Let us note that the 
(I) contribution contains single Sudakov logarithms, and in the (II) contribution (in i?7-boxes) there is also double 
Sudakov logarithms (DSL) but they are mutually cancelled out and reduced to single ones (see Section IV). Finally, 
the (III) contribution contains the DSL, and as it will be shown below, they predominate in the region M ^> mz- 
Besides, all of the contributions (I-III) contain zero power of Sudakov logs, i.e. terms without weak boson masses. 

Obviously, the collinear logarithms of (II) QED radiative corrections can compete with the DSL in the investigated 
region. This important issue has been studied at the one-loop level in || , where both the QED and weak corrections 
have been calculated for M < 2 TeV, but it has yet remained unsolved in the region of M > 2 TeV. Other important 
contributions in the investigated reaction at high invariant masses are the higher-order corrections (two-loop elec- 
troweak logarithms, at least), which also have been studied in the works El (see also the numerous papers cited 
therein) . Weak boson emission contribution has been recently calculated in [lfj and the contribution of higher-order 
corrections due to multiple photon emission has been computed in [TJ , these contributions are beyond the presented 
calculations. 

As has already been mentioned, for the future experiments at LHC aimed at the searches of NP in the reaction 
(fT|) it is important to know exactly the SM predictions, including the radiative background, i.e. the processes, which 
are experimentally indistinguishable from (fTJ). The important task is the insertion of this background into the CMS 
Monte Carlo generators, which should be both accurate and fast. For them to be fast it is necessary to have a set of 
compact formulas for the WRC which have been obtained here. Though we use the Asymptotic Approach (AA) (the 
Sudakov logarithms as the parameters of expansion in the investigated region of large M) , retaining the first and zero 
powers in the expansion we get good accuracy of calculation. 

Thus, we further introduce basic notations (Section II), show how to calculate the cross section corresponding to 
the (III) and (I) contribution diagrams (sec Fig.l,b-h): the heavy vertices (Section III), the heavy boxes (Section IV), 
and the boson self energies (Section V), we compare our results with those available at the quark level, and we give 
the numerical estimations of weak radiative corrections to the Drell-Yan cross section in the region 0.5 TeV < M < 
10 < TeV (Section V). 

II. NOTATIONS AND THE BORN CROSS SECTION OF THE DRELL-YAN PROCESS 

The Born cross section for the inclusive hadronic reaction AB — > l + l~ X is given in the quark parton model by 
formula 

a (M, yi () = l^-Re £ D"Dt>* ]T (t* + X u 2 )\^ £ F«{x + ,x.)\^ . (2) 

i,j="f,Z X— +- — q=u,d,s,... 
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Our notations are the following (see Fig.l,a): p\ is the 4-momentum of the first unpolarized (anti)quark with the 
flavor q and mass m g ; p 2 is the 4-momentum of the second (anti)quark with the flavor q and mass m q ; ki(k 2 ) is the 
4-momentum of the final charged lepton with the mass m; q — k\ + k 2 is the 4-momentum of the z-boson with 

the mass rrii\ Pa(b) is the 4-momentum of initial nucleons A(B). The invariant mass of dilepton is M — \[q 1 ~ . We 
use the standard set of Mandelstam invariants for the partonic elastic scattering s, t, u: 

s = (pi+ p 2 ) 2 , t = (p x - fa) 2 , u = (Ai - P2 ) 2 (3) 

and the invariant S — (Pa + Pb) 2 for hadron scattering. The first summation runs over two possible intermediate 
bosons: 7, Z. The third summation runs over all contributing parton configurations. The number 1/3 is a color factor, 
a is the electromagnetic fine structure constant. The propagator for j-boson has the form 

D]S = o 1 . =- , (4) 

s — rrij + irrijL j 

where Tj is the j-boson width. The combinations of parton density functions look like 

F q ± ( Xl ,x 2 ) = ff(x 1 )f?(x 2 )±ff(x 1 )f?(x 2 ), (5) 

where (x) is the probability of constituent q with the fraction x of the hadron's momentum in hadron H finding. 
The combinations of coupling constants for /-fermion with i- (or j-) boson have the form 

A//' = v)v) + a)a J f , X f y = v)a) + a)v f , (6) 

where 



v 



If ~ 2s 2 w Q f _ z If 



f ' f ' f 2swcw ' 2swcw ' 

Qf is the charge of fermion /, 7§ is the third component of the weak isospin of fermion /, and sw (cw) is the 

sine(cosine) of the weak mixing angle: sw = yl - c w , cw = mw /raz- 

Reducing the phase space of reaction (fTJ) to the variables M and y - rapidity of dilepton (we need also the variable 
t which is determined by r 2 = q 2 /S) we get a three-, two- and one-fold cross sections (without any experimental 
restrictions yet) 



dMdy 

and 



, 2 + 1 

a(M,y) ee — °— = / dC a(M,y,(), (9) 



-l°g# +1 



a(M) ee ^ = J dyjdC a(M, y, Q, (10) 

here the invariants s, i, and u mean 

s = M\ t = - l -s{\ - C), u = -~s(l + 0, (11) 

£ is cosine of the angle 9 between p\ and k\ in the center mass system of hadrons (C = cos 9) and the arguments of 
parton distribution functions in © have the form 

x± = re ±y . (12) 
Integrating over the whole region of £ we get much simpler expression 

a (M,y) = ^^-Re ]T D is D^*X^ ]T F*(x+, x-) A,*/. (13) 

i.j—j.Z q—u,d,s,. . . 
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III. HEAVY VERTICES 

To construct the contribution of Heavy Vertices (HV) to the radiative corrections (see diagrams (b-e) in Fig.l), we 
used the t'Hooft-Feynman gauge and an on-mass renormalization scheme which uses a,mw i^z-, Higgs boson mass 
win, and the fermion masses as i ndep endent parameters. Appropriate (for ultrarelativistic limit we are interested 
in) results can be taken from Ref.[12j|. where renormalized gauge boson fermion vertices for on shell fermions have 
been obtained. These results are presented as the form factor set to the Born vertices, so we can easily use them to 
construct the cross section: hence we replace the coupling constants in the Born vertex for the corresponding form 
factors: 

vj -^Fjf,a j f ~^F{ f . (14) 
Then the cross section of heavy vertices contribution will look like 

av(M,y,() = ^ Re ]T D is &** £ {t 2 + X u 2 ) ]T F«(x + , X-){\ F q ^\% + \f (15) 

"y,Z X— + i — q—u,d 7 s,... 

where the combinations of couplings constants and form factors are 

XJl j = F^v) + r/n' r \f* = 4/0} + F%v j f , (/ = q, I). (16) 

Electroweak form factors Fy* in ultrarelativistic limit have the form (there is the perfect coincidence with the 
formulas of Appendix C of Ref . [13| ) : 



F$ = ^[((vf) 2 + (af) 2 )A 2 (m z ) + ^-A 3 (m w )], (17) 



F} 1 = ^[2vfafA 2 (m z ) + -^A 3 (m w )], 
4tt As w 



F ^ = ^[((^) 2 + («^) 2 )A 2 K) - rVA 2 (m w ) + -^A 3 K)], 
4tt 2s^r 4s^ 



1 9 
-5-A 2 (m w ) + — J- 



F T = -^-[2^a|A 2 (mz)-^ 5 -A 2 (m w ) + ^2-A 3 (m w )], 



F v U = +^)')A 3 (m z )--^ r A i (m w ) + ^ r A a (m w )] ) 



t ,,,i 1 q 

FT = -^[2v z a z u A 2 {m z )-^ r A 2 {m w ) + —A 3 ( mw )]. 



F? = ^[vf((vf) 2 +3(aff)A 2 (m z ) + -^—A 2 (m w )-^A 3 ( mw )l (18) 
47T Ss^cw 4s^ 



^ = f [af(3(vff + (aff)A 2 (m z ) + -J—A 2 (m w )-^ 



F$ d = f [v z d {{v z d f + 3(a z f)A 2 (m z ) + 1 f u ° w A 2 (m w ) - ^K)], 
4vr Ss-^cw 4s^ 



F z A d = f [of (3(t>f )' + (aff)A 2 K) + 1 A 2 (m w ) - ^A 3 (m w )}, 

An as^cw 4s^ 



F^ = f [v Z ((v Z f + Ha z J)A 2 (m z ) - 1 ± \ Q d * w A 2 (m w ) + ^A 3 (m w )}, 
47T os^cw 4s^ 



FT = f [^flfrf )' + (a z S)A 2 ( mz ) - 1 ± f dS " w A 2 (m w ) + |g^A 3 (m w )]. 
4vr Ss^Ch/ 4s^ 
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Here it is necessary to note that in on-mass renormalization scheme the self-energies of it-quarks' diagrams also 
give a non-zero contribution to the cross section. However, the results for it-self-energy are factorized in the same 
manner as the vertices. It gives us a possibility to sum both contributions in one formula, where the u-self-energy 
contribution is completely cancelled with the corresponding terms of heavy vertices. So, it could be said that the 
formula (|15[) gives the cross section induced by both the heavy vertices and self-energy of it-quarks. 

Exact expressions for A2,3 can be found in [l2l ]. and there we give the real part of expressions for functions A2,3(mi) 
that provide a good accuracy in large dilepton mass region (see Introduction): 

7T 2 7 5 1 

A 2 (m 4 ) = — - - - 3Z i;S - Z? s , A 3 (mi) = - - -Z i)S . (19) 
The Sudakov logarithm that we denote as U yX has the form 

o 
TTl • 

k,* = log t-t (i=Z,W; x = s,t,u). (20) 



IV. HEAVY BOXES 



The calculation of a two heavy boson (box) contribution is more complicated procedure since it requires the 
integration of 4-point functions with the complex masses in unlimited from above kinematical region of invariants. 
The difficulties of such a procedure have been pointed out as far back as a pioneer paper, Ref. Fortunately there 
is a way to avoid many troubles with the integration all of the terms in the box contribution using a rather simple 
algebra. Let us explain this way. 

First of all we construct the ZZ-box cross section for qq — > (see diagrams in Fig.l,f-g) using the standard 

Feynman rules: 



its q (2tt) 
here the 2-particle phase space element is read 



da zz = ^dT.Rej^ I d*k £ D*°*(D™ + C zz ), (21) 

fc=7,Z 



d 3 ki d 3 k 2 

and D ZZ (C ZZ ) is the direct (crossed) diagram contribution. 

Neglecting the fermion masses we present the direct contribution in the form: 

D zz = TI&TrWh'Yu&i ~ k) 7v pf' k ( Pl )]Tr[ loi k 2 ^(k - *i)7"pf Z,k (*i)], (22) 
where the integrand of 4-point scalar function looks like 

r\D \ (OQl 

ZZ ((q- k) 2 -m 2 z )(k 2 -m 2 z ){k 2 ^ 2k l k){k 2 - 2p x k)' ' 

Combinations of the density matrices p(p) and the coupling constants can be reduced to the production of A-factors 
(see formulas ©) and p = J^p^ (here p — pi,P2, k\, ki or k) 

pf *(P) = («f - af l 5 )p(p)(v k f + af l5 ) = h\ f f* - X f zz ^ 5 )p, (24) 



v 



f = (vf)+(af), af = 2vfa z f . (25) 



To extract the part of cross section which predominates in the region \x\ 3> m 2 z (see x in formula (I20p ~) we should 
make the equivalent transformation of the cross section based on the close connection of infrared divergency cross 
section terms and the Sudakov leading-log terms: 



(26) 
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The integral over k of the first term in (|26[) is 



(2tt) 2 7 " q 2 -m\" (2ir) 2 J k 2 - m% L(fc 2 - 2k 1 k)(k 2 - 2 Pl k) ' (k 2 - 2k 2 k)(k 2 - 2p 2 k) 

here we used the trivial correlation typical for the Born cross section 

S ZZ ' k = \^b a PWf ^(piWlak^pf^ih)] « bf'H 2 + fc^V, (28) 

where 

Then we integrate one of terms in (|27p over fc using the standard method of paper [l4| (another term is integrated 
similarly) : 

_l} j dy 



(2tt) 2 J (k 2 - m 2 z ){k 2 - 2fc 2 fc)(fc 2 - 2p 2 k) 4 J J m%(x - y) + (p 2 (l - x) + k 2 yf 



o o 

i 



1 f, 1 , t(x-l) 1 ,tt 2 1 , 2 -t - mi T 

- / da;- r 5- log- — 5—^ = - — ( — + - log 2 g— ^ + Li 2 , z ), 30 

A J t(x-l)-m% m| 4i v 3 2 5 m| m| + r 



here Li2 denotes the Spence dilogarithm. Retaining the terms which are proportional to the second (~ l 2 x ), first 
(~ l\ x ) and zero (~ 1® X = 1) power of the Sudakov logarithms and thereby neglecting the terms which are inessential 
in the region of large dilepton mass \x\ m 2 z we get the asymptotic expression 

d A kDl z ^--B' LZX {^ + l 2 zt ). (31) 



(2tt) 2 J 1 s v 3 

Surely in these (and subsequent) expressions we can retain only leading ~ if x term as it has been done, for example, 
in [l5j . Naturally in such an approximation our results coincide with the results of that paper. However we are able 
to retain the first and zero powers of expansion leading-log parameter and we will use this opportunity to improve the 
accuracy of radiative correction estimation (one remark more - we retain the lj x and 1® X also in the heavy vertices 
part). 

The asymptotic expression (at |x| 3> m z ) for the second term of the formula ([26]) can be presented as 

DT « ^ ZZ ln^n ?T - 4n£ (g -^ 2 + *% B zz * (32) 
Reducing the vectorial and tensor integrals to the scalar ones we get 

d^kDf » -2fe zz ' k [(Gjr + G^)(t 2 + u 2 ) - 2(R - N )u - X Q t^^-} - ibf' k t 2 [G^ + - X -](33) 



(2tt) 2 J z ~ -~ LV ~ U ' ~ u "~ '~> ■■">- s ' ™+ ~ L ~ u ' ~ u "V 

where scalar integrals are 

C m = -*— f - 1 Ho- 2 q2 k— 1 

(2tt) 2 J (k 2 -2k 1 k)k 2 (k- q) 2 8q 2K & m 2 3 h 

,M i f d 4 k 1 fl _ 2 g 2 



'~ J (k 2 -2 Pl k)k 2 {k-q) 2 ^^Zq 2 [ ° g m| T j ' 

i f d 4 k 1 g 2 

* = P(fc-,) 2 = 4 (1 ° g L 

i /" d A k 1. m 2 m |2fcipi| 

= (^py (fc 2 ~2k 1 k)(k 2 -2 Pl k) = 4 ( g X~ g mT + g "mmT^ 
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(2tt) 2 7 (fc 2 - 2k 1 k){k 2 - 2 Pl k)k 2 (k- qf 



1 1 |2pifci|m 2 . |2piA;i| , 1. \2pih\m 2 \2 Pl h\ tt 2 



v / (o lo s 3 — lo s — r~ + o log 3 — log — 2 5-)- ( 34 ) 

These integrals have been calculated on the analogy with the corresponding ones in paper of J.Kahane [16| (with 
the important distinctions due to different reaction channel: i-channel in 16j, and s-channel in this paper). The 
expressions (1341) contain the parameters of two sorts: masses of fermions and L - parameter regulating the ultraviolet 
divergence, both of them are completely cancelled out in the final expression 

• / d^Df « b zz > k l s u + (b zz > k (t 2 + u 2 ) + 2bf 'V)i/ S 2 , l s = log ±. (35) 
J 2s \t\ 



The contribution of crossed part (Fig.l,g) to the cross section can be calculated on the analogy with the direct one. 
Besides, there is an interesting symmetry between the direct and crossed parts (see, for example, [5j), namely 

C zz = ~D zz \ t ^ u ^ & zz,k^ b zz, k . (36) 

And now let us present the total final expression for the contribution of ZZ-box into cross section of the Drell-Yan 
process at large invariant dilepton mass: 

2a 3 



+f£{x+)f* (x-)(S ZZ ' k (u, t, &+,&_) - S zz > k (t, u, b-,b+))] , (37) 



3SM , 

h— 7,z q—u,a,s, 



where 



oZZ.k r>_2 ;2 

«5 zz ' k (i, u, b+, b-) = + l%, t ) ~ b zz - k ul s - (b ZZ ' k (t 2 + u 2 ) + 2b ZZM t 2 )±. (38) 



Integration over £ of the formula (|38[) gives 



i +i 

ZZ,k/ + „. i. . \ j/- _ / cZZ.k/ 



<(t,u,b+,b-) d( = J 5 LL ' k {u,t,b + ,b-) d( = 
-l -l 

= £(l2A4 z ' k A ; zz ' k (^ + l\ s ) + (226 zz < k + 46 zz ' k )^ + 276 zz ' k ) (39) 

(it is convenient to use it for the construction of twofold cross sections crzz,ww(-^j y))- 

To obtain the WW-box contribution into the Drell-Yan cross section using the expressions (|37|) and ([38| one 
should: 1) do the trivial substitution Z — + W in all indices of coupling constants and boson masses, 2) take into 
consideration that some quark diagrams with two W-bosons are forbidden by the charge conservation law. Then 
WW-box contribution will look like 



aww(M,y,C) = ^Re £ D ks * (x_)£ ww ^, u, 6 +1 6_) + /^(x + )/f(^)<5™( U , t, b+, &_)' 

k=~y,Z q=u,c 

- E [ft^+)f q 3 ^-)S WW '\^t,b-,b + ) + f^(x + )ff(x^5 ww ' k (t,u,b^b + )]). (40) 

This second feature of WW-boxes explains the domination of WW-contribution into the Drell-Yan cross section in 
comparison to ZZ-contribution (see the section on numerical analysis). The point is that the ZZ-box cross section is 
proportional to the difference 

5 zz - k (t, u, b+, 6_) - 6 zz - k (u, t, 6_, 6+) ~ l\ t - l\ u = log|(& )t + J, (41) 
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(see also [15|]) i.e. it ~ lzx> whereas the leading parts of WW-cross section do not contain the difference (|4"Tj) and 
are proportional to 1^ . Let us note here that the factorization property (|41[) is absent in heavy vertex part and is 
present in infrared finite part of "/Z-box contribution, which can be easily obtained on analogy with ZZ-box using 

<F z ' fc (t,it,&+,6_) = + l%,t) ~ 2bl Z ' k ul s - (bl z < k (t 2 +u 2 ) +2bf' k t 2 ) l X (42) 

We attribute the 7Z-DOX to the "other QED corrections" (see (II) contribution in Introduction) and do not consider 
it here any more. 

V. DISCUSSION OF NUMERICAL RESULTS 
A. Numerical estimation at the parton level and comparison with the available results 

Now we discuss the observable quantities at the parton level. The reason of it is twofold: to study the characteristics 
of the WRC in the absence of parton distribution functions and (what is much more important) to compare them 
with the available results. Except our paper [Tt]], where the WRC to Drell-Yan production have been calculated by 
different methods: exact integration, numerical integration over Feynman parameters [Lsl ] and the AA considered 
here, there is the pioneer paper Q and the calculation of SANC group [ill. 

At the parton level there is a possibility to compare the results with Q, since in this paper the set of relative 
non-QED corrections has been presented at high parton center of mass energies for the total cross section (FIG. 3 of 
H) and for the parton forward-backward asymmetry (FIG. 4 of 5]). Further we present our variant of corresponding 
WRC in the [5J region (up to 1 TeV) and in the region of higher energies (up to 10 TeV). For the correct and tuned 
comparison we used in this subsection the same set of the SM parameters which was used in [f| and all the rest 
prescriptions of this paper. We also add the contribution of boson self energies to the numerical estimation (see 
Fig.l,h) 

i,j=7,Z X=+>- 

+w z D Zs EH-* (W/^ + X «x j h™){t 2 + X u% (43) 

i=J,Z X = + ,~ 

Here n 7,Z ' 7Z are connected with the renormalized photon-, Z- and 7Z-self energies [12J as 

V7 yZ y~fZ 

if = — , n z = if z = — . 

s s — m z s 

Thus, Fig. 2 shows the relative corrections to the total cross sections (the left picture) and forward-backward 
asymmetry at the parton level (the right picture) for uu — > /i + fi~ as the functions of \fs. As in [f|, for e + e~ final 
state we obtain the identical results. The distinction between the dd- and ira-scattering is the same as in [f|, therefore 
we do not present the graphs for dd-case here. The symbols BSE, HV, ZZ, WW in Fig. 2 denote the boson self 
energy, heavy vertex, ZZ-box, WW-box contributions correspondingly, the dashed line corresponds to the sum of all 
contributions. The asterisks are the points taken from FIG. 3 and 4 of paper [5(. 

A rather good coincidence with the presented results and the results available from || takes place: we can see 
in Fig. 2 almost the same scale and behavior of relative corrections (the decrease with energy, the convexity and 
t-quark threshold kink) and the tendency of the results to converge with the increase of energy. We remind that 
the approaches in both calculations are not quite similar. Although for the calculation of the BSE and HV parts 
here and in Ref.[5J the same BSE insertions and form factors have been used, but the heavy box contribution has 
been calculated in different ways: via the direct 4-point functions integration in [5J, and using asymptotic property 
of heavy boxes expressed through the Sudakov logarithms powers in this paper. It is obvious that the difference of 
the corrections in different approaches strongly depends on the scale of energy: the asymptotic Sudakov expansion 
works well at large values of y/s only. Fig. 2 shows that at yfs=l TeV this difference is rather small both for total uu 
cross section (~ 1.5%) and for quark forward-backward asymmetry (~ 0.001) and with the increase of yfs is found 
to become still smaller. It is a pity, but in 5] the parton observables in the region yfs > 1 TeV have not been 
investigated. 

To compare_our results with the available in the region of high invariant masses y/s > 1 TeV we contacted to 
SANC group [HI and to the authors of 0] (see also program ZGRAD [1(|). Having used the same sets of electroweak 
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FIG. 2: The relative corrections to the total cross section (the left picture) and the forward-backward asymmetry at the parton 
level (the right picture) for uu — > /i + fi~ as a functions of yfs. The symbols BSE, HV, ZZ, WW denote the boson self energy, 
heavy vertex, ZZ-box, WW-box contributions correspondingly. The dashed line correspons to the sum of all contributions. The 
asterisks are the points taken from FIG. 3 and FIG. 4 of paper 

parameters we got excellent coincidence for all WRC (see Table 1), which shows the same (as in Fig. 2) corrections 
(in per cents) to im-cross section as a functions of yfs in the region 0.1 TcV < yfs < 10 TeV. 
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-3.107 


-4.8739 


-1.6883 


-1.6688 


2.4481 


11.2259 


12.2144 


0.5 


-0.2144 


-0.2145 


-0.1895 


-10.777 


-10.1087 


4.2958 


4.2978 


4.0943 


11.1526 


11.9455 


1.0 


-0.3346 


-0.3346 


-0.3251 


-16.998 


-16.5720 


6.2447 


6.2451 


5.9910 


12.2096 


12.9793 


2.0 


-0.4638 


-0.465 


-0.4612 


-25.442 


-25.2497 


8.5534 


8.5535 


8.4247 


13.1993 


13.9634 


3.0 


-0.5423 


-0.543 


-0.5410 


-31.468 


-31.3554 


10.1709 


10.1710 


10.0935 


13.7682 


14.5314 


5.0 


-0.6432 


-0.643 


-0.6416 


-40.185 


-40.1306 


12.4894 


12.4895 


12.4512 


14.4811 


15.2437 


10.0 


-0.7787 


-0.779 


-0.7782 


-53.989 


-53.9695 


16.1160 


16.1160 


16.1024 


15.4456 


16.2080 



Table 1. The relative corrections (in percents) to the ZZ, WW, HV and BSE cross sections at the parton 
level for uu — > as functions of yfs, calculated by different groups: SANC[19|, program ZGRAD [20j and using 

the results of this paper. 



Let us analyse the ZZ-box: we can see that in the region of not so large values of yfs (0.1 - 0.2 TeV) the results 
are rather different (although the corrections hold the similar behaviour and scale), but starting with yfs ~ 0.5 TeV 
the difference is small and tends to zero with the growth of yfs. Of course, this fact confirms the precision of both 
calculations and estimates the region where the asymptotic approach suggested here works well. I suppose this point 
is 0.5 TeV, here the relative error of asymptotic (to exact SANC) calculation is ~ 0.1161, whereas at 1 TeV (2 TeV) 
the error is ~ 0.0284 (~ 0.0056). 

There is no special need to compare separately the direct and crossed parts of ZZ-box correction (the same way 
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as WW-box) containing DSL, since having the good agreement for the whole ZZ-box, which is proportional merely 
to the first power of Sudakov logarithms (see we will obtain the agreement of the same order. Nevertheless, 

we compare the results for the direct WW-box: for the point 0.5 TeV the relative error of asymptotic (to SANC) 
calculation is ~ 0.0620, and at 1 TeV (2 TeV) this error is ~ 0.0251 (~ 0.0076). Table 1 also shows good agreement 
for HV part of radiative correction: the first HV column is the result of SANC, the second one is obtained using the 
formulas of Section III and exact expessions for A2,3 from [121 ] . the third column corresponds to the approximation 
(I19[) : we can see a very good coincidence of the first and second columns and a rather good agreement between the 
first and third columns - the relative error is 0.040 at y/s=l TeV and 8.44 x 10~ 4 at ^=10 TeV. For the BSE part 
the agreement is not so good as for HV - there is a small constant difference ~ 0.5% obviously caused by the use of 
various gauges. 

Let us pass to the numerical analysis of hadron observables. 



B. Numerical estimation for future CMS program 



Futher the scale of weak radiative corrections and their effect on the observables of the Drell-Yan processes for 
future CMS experiments will be discussed. We used the following set of electroweak parameters (such as in |20jj): 
a = 1/137.03599911, m w = 80.37399 GeV, m z = 91.1876 GeV, the energy of LHC y/s = 14 TeV and the CTEQ6 
set of unpolarized parton distribution functions (PDF) [2l[ (with the choice Q 2 = M 2 ). 

We impose the experimental restriction conditions on the detected lepton angle —(* < C < C* an d on the rapidity 
|j/(0l — 2/(0*- For CMS detector the cut values of (* and y(Z)* are determined as 

0* 

y(l)* = - log tan — = 2.4, C* = cos 6* w 0.9837. (44) 

Then the expression (|10[) with consideration of the experimental restrictions (|44p and the use of theta-functions can 
be modified to 

cr(M)= J dy J d( a(M, y, Q9(cosa + C*)#(— cos a + (*), (45) 

where a is the scattering angle of the lepton with the 4-momenta fc 2 (a = Pi, k 2 ) in the center mass system of hadrons. 
This angle has a non-trivial relation with 9 and y: 

cosfl-/ . /sing e 2 y - 1 

a = 7r — arccos — ■ — arcsm — , , f = — . ( 46 ) 

^1 + J 2 -2/ cos (9 y/l + f 2 -2fcos6 e 2y + 1 

Let us note here that the second standard CMS restriction pr{l) > 20 GeV in the region of large M considered here 
and at (|44[) is satisfied completely and automatically. It can be seen in Fig. 3 where all of the CMS experimental 
restrictions for M = 1 TeV have been shown. 
In Fig. 4 the relative corrections are presented 

c _ a c (M) 

hl ~ a (M) ' (47) 

as functions of M in the region 0.5 TeV < M < 10 TeV with experimental restrictions of the CMS. Index "C" in 
all numerical analysis means 

C = ZZ,WW,HV,BSE,tot; tot = ZZ + WW + HV + BSE. (48) 

It can be seen, that the BSE, HV, and WW-contributions are most significant, and the WW boxes give the dominant 
(negative) contribution. In Fig. 4 we can also see the significant contribution of zero power of Sudakov logs (~ lwx) : 
the difference "WW (DSL only)" - "WW (total)" is ~ 4% and depends weakly on M. Big contribution of lyy x justifies 
the detailed calculations which have been done in Section IV. The 4%-difference is caused mainly by the 27r 2 /3-term 
in (31]), the contribution of Df w (32]) is far less than 1%. The total WRC effect is significant: in region M < 1.5 TeV 
is positive, since the sum of positive HV and BSE contributions exceeds the negative WW-boxes, but for M > 2 TeV 
the total correction is definitely negative and decreases with the increase of M. Thus, the S^f correction changes the 
dilepton mass distribution up to ~ +3(— 12)% at M = 1(5) TeV. 
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FIG. 3: The region of the cross section integration over £ and y taking into consideration the CMS experimental cuts at 
s/S = 14 TeV and M=\ TeV. 




FIG. 4: The relative corrections 5 M corresponding to the ZZ, WW, HV, BSE-contributions with experimental restrictions of 
the CMS as a functions of M. The dashed line corresponds to the sum of all the WRC contributions. 



Here we should note that the importance of the investigated problem demands the cross-checking of this part with 
the earlier obtained in [|| . The direct comparison between numerical results of the two works is not possible, since in 
Ref.Q the effective Born approximation has been used to calculate the relative corrections. We do not use the effective 
Born approximation as the main subject of our investigation is the region of large invariant mass (M > 1 TeV) where 
this approach is illegal. Fortunatly, recent paper [l(| provides the possibility to compare the results of Ref. [j| and 
presented by us in the region of high invariant dilepton mass (see the first picture in Fig. 7 of [HI). At first glance 
the results look absolutely different: for M = 1(2) TeV the relative corrections of is ~ — 5%(— 11%) whereas our 
result is ~ +3%(— 2%). However, the solid line in Fig. 7 of correspons to total one-loop result and to compare 
correctly we should add the photon corrections: "soft" and "hard" photon contributions. We have our own code for 
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calculation of this part (FORTRAN program READY - "Radiative corrEctions to lArge invariant mass Drell-Yan 
process" , which is supposed to be published in the near future, for more details see [22| ). Main features of READY are 
the same as in Ref.[3]: using the " soft" -" hard" separator u and independence of total result on it, lepton identification 
cuts reduce extremely the "hard" FSR contribution and so on. For the decision of quark mass singularity problem 
we used the procedure of linearization [19(. The "soft" part, which we should add to pure WRC for comparison 
with [io| . is large and negative (for example, at M = 1 TeV, u> = 10 GeV and I = e the FSR-part of "soft" relative 
correction defined by formulas (72), (73) in Ref-fl^ equals —40.7%. The positive "hard" FSR-part suppressed by 
lepton identification requirements will be far less, and the rest parts (ISR after quark mass removing, and the mass- 
independing interference) will be small as compared to the FSR. Hence, we are supposed to have a rather good (at 
least qualitative) agreement with hadron estimations of Ref. [l(| : the relative correction will be negative and decrease 
with M in region M > 0.5 TeV. Unfortunatly the READY is not quite complete and here we can not present more 
exact numbers, but we hope future publications of the results of code READY must clarify this problem later. 

In the end of this section we want to discuss the effect of the PDFs on the weak radiative corrections in this very high 
mass regime, and what uncertainty they introduce. Comparing the results for CTEQ6 set [2l[ and MRST 2004-QED 
set [23| of PDFs the essential effect of PDF choice for the cross sections has not been noticed, the relative error is on 
one per cent level and has no tendency to increase with the increase of M, at least at M ~ 3.16 TeV (unfortunately 
MRST 2004-QED set for Q 2 > 10 7 GeV, i.e. for M > 3.16 TeV, does not work, and it does not allow to compare the 
result for both sets of PDF in the region of extremely high invariant masses). 



VI. CONCLUSIONS 



The weak radiative corrections with the neutral current to the Drell-Yan process for large invariant mass of a 
dilepton pair have been studied. The compact asymptotic expressions have been obtained, which expand in the 
powers of the Sudakov electroweak logarithms. At the parton level the investigated radiative corrections have been 
compared with the existing results and a rather good coincidence at energy > 0.5 TeV have been found. The numerical 
analysis in the high energy region has been performed with use of the standard CMS cuts. It has been found that 
the considered radiative corrections become very significant at high dilepton mass, and their large scale does not 
permit to neglect this part of the radiative correction procedure in the future experiments at LHC directed on the 
NP research in the Drell-Yan process. Some part of corrections has become beyond the scope of the presented paper 
(the QCD corrections, the pure QED corrections, and the two-loop electroweak logarithms), and still it requires a 
thorough analysis in order for us to completely solve the radiative corrections problem to the Drell-Yan process at 
extremely high energies. 
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